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We predict the enhanced transmissivity of modulated slabs of layered superconductors for ter- 
ahertz radiation due to the diffraction of the incident wave and the resonance excitation of the 
eigenmodes. The electromagnetic field is transferred from the irradiated side of a slab of layered 
superconductor to the other one by excited waveguide modes (WGMs) which do not decay deep into 
the slab, contrary to metals, where the enhanced light transmission is caused by the excitation of the 
evanescent surface waves. We show that a series of resonance peaks (with T ~ 1) can be observed in 
the dependence of the transmittance T on the varying incidence angle 0, when the dispersion curve 
of the diffracted wave crosses successive dispersion curves for the WGMs. 
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Since the first observation by Ebbesen et al. [T], the 
enhanced light transmission (ELT) through metal films 
perforated by subwavelength holes is the focus of atten- 
tion for many research groups (see, e.g., reviews [3] 
and references therein). This phenomenon is observed in 
films with thicknesses much larger than the skin depth, 
and the transmission coefficient turns out to be much 
larger than that predicted by Bethe's theory of electro- 
magnetic diffraction at small apertures [4|. The ELT is 
related to coupling of surface plasmons resonantly ex- 
cited at both sides of the perforated film. Discussion of 
this and some alternative mechanisms of extraordinary 
transmission can be found in the review by Zayats et 
al [5]. Recent interest to the aforementioned effect is 
due to its possible applications for light control, photo- 
voltaics, detection and filtering of radiation in visible and 
far- infrared frequency ranges. 

For observation of the ELT in thick enough metal films, 
there should obviously exist a mechanism for the transfer 
of the electromagnetic energy from the irradiated side of 
the film to the other. In regards, we would like to men- 
tion two forgotten ways to make metal films transparent 
for electromagnetic waves. Both of them are character- 
istic for pure metals at low temperature. The first one 
is the so called anomalous penetration of the electromag- 
netic field deep into the film by the chains of the Larmor 
electron orbits (see, e.g., Ref. [B]). In the external dc 
magnetic field parallel to the sample surface, the elec- 
trons with long enough mean free path carry the elec- 
tromagnetic field out from the skin layer and form an 
additional current layer at the distance of the Larmor di- 
ameter from the sample surface. The another group of 
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Larmor electrons carry the electromagnetic field deeper 
into the sample and form the next current layer, and so 
on. Thus, the Larmor electrons can serve as carriers for 
the transfer of the electromagnetic energy through the 
metal films. 

The second mechanism is related to different kinds 
of weakly decaying electromagnetic waves: helicons, 
dopplerons, cyclotron waves, etc. (see, e.g. Refs. [IHS]). 
Under definite special conditions, in the presence of the 
external dc magnetic field, they can propagate in a metal 
and carry the electromagnetic energy from the irradiated 
side of a film to the other. 

In this Letter, we consider a novel mechanism for the 
ELT, which is a specific combination of the two mecha- 
nisms discussed above. We predict and analytically study 
the enhanced transmission of terahertz radiation through 
slabs of layered superconductors with periodically mod- 
ulated values of the maximum density J c of the c-axis 
Josephson current. Such a modulation can be realized ei- 
ther by irradiating a standard Bi2Sr2CaCu20g+,5 sample 
covered by a modulated mask |10j or by pancake vortices 
controlled by an out-of-plane magnetic field [TT]. Simi- 
larly to the case of the ELT in metals, the modulation 
results in the generation of diffracted waves, which, un- 
der the resonance conditions, excite the electromagnetic 
eigenmodes. 

The Josephson current along the crystallographic c- 
axis couples with the electromagnetic field inside the 
insulating dielectric layers, forming Josephson plasma 
waves (JPWs) (see, e.g., review [T2] and references 
therein). Thus, the propagation of electromagnetic waves 
through the layers is favored by the layered structure. 
The study of these waves is very important because of 
their terahertz frequency range, which is still hardly 
reachable for electronic and optical devices. It is nec- 
essary to emphasize a principal difference of the JPWs 
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in the strongly anisotropic Josephson plasma of lay- 
ered superconductors from the electromagnetic waves in 
isotropic media, e.g., in standard metals. Contrary to 
waves in metals, the dispersion law for JPWs has a "hy- 
perbolic" form [T2"] . 
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Here uj is the wave frequency; k s x and k s z are the compo- 
nents of the wave vector along and across the supercon- 
ducting layers; A c = c/wje 1 / 2 and A a b are the magnetic- 
field penetration depths along and across the layers, re- 
spectively; ujj r*> 1 THz is the Josephson plasma fre- 
quency; £ is the interlayer dielectric constant; c is the 
speed of light. 

Equation shows that JPWs can propagate across 
the layers for u> > ujj only. Moreover, the longitudi- 
nal component k sx of the wave vector should be large 
enough, 



k s x > k c = — {uj/uJj) 2 - 1. 
An 
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This specific feature of the JPWs is a reason for the 
principal difference of ELT in metals and in layered su- 
perconductors. For metals, not only the basic wave 
with k sx = (uj/c) sin 6, but all diffracted waves with 
k sx — (uj/c)sin9 + ng, exponentially decay to the mid- 
dle of a sample (here 9 is the incidence angle, g is 
the period of the reciprocal lattice of modulations, n 
is an integer). Contrarily, for layered superconductors, 
while the basic wave with (cj / c) sin 8 < k c exponen- 
tially decays into the sample, the diffracted waves hav- 
ing \(uj/c) sin 9 + ng\ > k c , can propogate across the lay- 
ers and, therefore, can serve as carriers for the electro- 
magnetic energy, similarly to the helicons, dopplerons, 
and cyclotron waves described in Refs. [7H9]. Thus, the 
diffraction in metals can result in the resonant excitation 
of the symmetric or antisymmetric evanescent surface 
waves, whereas the diffraction in the layered supercon- 
ductors can provide the resonant excitation of waveguide 
modes (WGMs) (see Refs. OH!]), which do not decay 
deep into the slab but oscillate across the layers. 

In this Letter, we analytically study the ELT in lay- 
ered superconductors for the simplest case of weak har- 
monic modulation under conditions when the inequalities 
(u/c) sin 9 < k c < \(oj/c) sin 9±g\ are fulfilled. Thetrans- 
mittance T versus the incidence angle 6 displays a series 
of resonance peaks with T ~ 1 emerging by the successive 
intersection of different dispersion curves for the WGMs 
with the dispersion curve of the diffracted wave. 

Electromagnetic field in vacuum and layered supercon- 
ductor. — Consider a slab of layered superconductor of 
thickness d surrounded by the vacuum (see Fig.[T]). The 
crystallographic ab-plane coincides with the a;j/-plane 
and the c-axis is directed along the z-axis. The plane 



z = is the top edge of the slab. Suppose that the max- 
imum c-axis Josephson current density, J c , and, there- 
fore, the Josephson plasma frequency ujj are periodically 
modulated in the x-direction with a spatial period L, 

loj(x) = ujj [1 + / cos( 5 a;)] , g = 2ir/L, f < 1. (3) 
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FIG. 1: (Color online) Schematic geometry of the problem; k ! , 
k r , and k* are, respectively, the wave vectors of the incident, 
reflected, and transmitted waves; qi is the longitudinal wave 
vector of the excited waveguide mode (WGM). 

A plane electromagnetic wave of TM (transverse mag- 
netic) polarization, E inc = {E x nc , 0, E™ c } and H inc = 
{0, -ff lnc ,0}, is incident from the vacuum onto the top 
side, z = 0, of the superconductor at an angle 8. The 
in-plane and out-of-plane components of its wave vector 
k l are 

k x = q = k sin 9, k z = k cos 9, k — ujje. (4) 

The in-plane periodic modulation results in generating 
the diffracted waves with the in-plane wave numbers q n — 
q + ng. For simplicity, we derive the electromagnetic 
field distribution taking into account only the first-order 
diffracted wave, q\ — q+g, which provides the ELT under 
conditions when q\ is close to the wave vector of one of 
the eigenmodes. Then, the results are generalized to the 
minus-first diffraction order, g_i = \q — g\. 

We assume that the diffracted wave in the vacuum 
is evanescent {q\ > oj/c), and its in-plane wave num- 
ber qi is close to the corresponding wave number of one 
of the eigenmodes. Thus, the electromagnetic field in 
the vacuum over the superconductor (z < 0) is pre- 
sented as a sum of the incident wave (with unit ampli- 
tude), specularly diffracted (n = 0) wave, and the first- 
order diffracted (n = 1) evanescent wave. The magnetic 
HY p(x,z) and the tangential electric E xtop (x, z) fields 
read 
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(x, z) — e~xjp(iqx + ikz cos ( 



i?o exp(iqx — ikz cos 9) + R\ expliqix + >cY z), 
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(x, z) = cos [exp(jgx + ikz cos ( 



i?o exptiqx — ikz cos 9)] ^— R\ expfzqix + xY z), 

k 
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where the attenuation coefficient x\ = \Jq\ — k 2 > 0. 

The electromagnetic field under the superconductor 
(z > d) is a sum of the transmitted wave and the first- 
order diffracted evanescent wave, 



H^ ot (x, z) — Tq exp[iqx + ik(z — d) cos 9] 
+ T\ exp[igix — x\{z — d)], 

E^ hot {x, z) = cos9T exp[iqx + ik(z — d) cos( 
■ V 

H — T\ exp[jgix — xY (z — d)]. 

k 



(7) 



(8) 



The electromagnetic field inside the layered supercon- 
ductor is determined by the distribution of the gauge- 
invariant phase difference tp(x,z, t) of the order parame- 
ter between the layers (see, e.g., Ref. [H]), 



dH s 
dx 
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Here, wj(x) is given by Eq. ([3]), ujj — (S-KeDJc/he) 1 / 2 
is the Josephson plasma frequency unperturbed by the 
modulation, J c is the maximum value of the Josephson 
current density j z = J c sin ip, D is the spatial period 
of the layered structure, $o = ffch/e is the magnetic 
flux quantum. The dimensionless relaxation frequencies 
r a b = 4ira a i,uiX 2 b / 'euj 2 X 2 and T c — 4ira c u! / 'euj 2 are pro- 
portional to the averaged quasiparticle conductivities a ab 
(along the layers) and a c (across the layers). Note that 
r c can actually be neglected due to the smallness of the 
out-of-plane conductivity a c . 

The phase difference ip is governed by a set of coupled 
sine-Gordon equations (see, e.g., review [12] and refer- 
ences therein). For linear JPWs in the continuum limit 
the coupled sine-Gordon equation can be written as 



1 - (1 + iT ab )X 2 ab 
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dx 2 
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= 0. 



Equation ( 10 1 can be solved perturbatively in small mod- 
ulation amplitude / < 1, see Eq. Then, Eq. |9} 
gives the magnetic field inside the layered superconduc- 
tor (0 < z < d), 



H s (x,z) = * (x) 

x [C£ exp[p (z - d)} + Cq exp(-po^)] 
+ *i(x) [Cf exp{ixf z) + Cf exp(-ixfz) 



(11) 



with 



*o(x) = exp(iqx) - F i exp(iqix) , 
*i(x) = exp(iqix) + F 01 exp(iqx) ; 
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The tangential component of the electric field in the lay- 
ered superconductor reads 

E x {x,z)=i(l + ^f) 

x |a * (x) [C$ exp[p (z - d)] - cxp[-p z]] ( u ) 



a% ^i(x) [Ci exp[i>rf z\ — C x exp[— ixf z]\ j 



where 



a = kX a b\ 1 - 
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are the small surface impedances for the basic and the 
first-order diffracted waves, respectively. 

Transmittance and reflectance. — Matching the tan- 
gential components of the electric and magnetic fields at 
the boundaries, z = and z = d, of the superconduc- 
tor, we obtain eight linear algebraic equations for eight 
unknown amplitudes, Rq, Ri, To, T%, Cq , , C± , and 
Ci . Solving these equations yields the transmittance 
| To | 2 and reflectance |i?o| 2 of the superconducting slab, 



|Th| = 



4 a\ I cos 2 9 
V 2 +B 2 ' 

V 2 + (xfdr o6 /2)' 
T> 2 + B 2 



where 



V = tan (>rf d) — 2 — y a\ 



B 



xfd 



2F 2 



ai 
cos 9 



(16) 
(17) 

(18) 
(19) 



Equation V = with / = 0, i.e., with F = in Eq. (13 1, 



defines the spectrum of symmetric and antisymmetric 
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waveguide eigenmodes [T3] . The weak modulation of the 
plasma frequency ^ results in a shift of the dispersion 
curves and in additional damping, namely, in the leakage 
of the eigenmode energy due to diffraction. This leakage 
is specified by the second term in Eq. (19). 
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Equations (16) and (17) describe the resonance en- 



hanced terahertz-wave transmissivity and, respectively, 
the suppression of the reflectivity due to the excitation 
of the waveguide mode by the first-order diffracted wave. 
Obviously, similar expressions can be derived for the 
case of the resonance waveguide mode excitation in the 
minus-first diffraction order. To do this, one should sub- 
stitute the wave number q_i — \ksui9 — g\ instead of 
qi = g + k sin 9 in all formulas. Figure [2] shows six reso- 
nance peaks in the dependence of transmittance |To| 2 and 
reflectance |-Ro| 2 on the incidence angle 9. Three of them 
are due to the first-order resonance diffraction when the 
argument xfd of the tangent in Eq. (18) is close to mn 
with m = 8, 9, 10; and the others are related to the reso- 
nances in the minus-first diffraction order at x^d ~ rmr 
with m — 4, 5, 6. 

We also illustrate the ELT effect by the distribution of 
the magnetic field in Fig. [3j Under the nonresonant con- 
ditions, the interference pattern is seen in the vacuum re- 
gion over the superconductor, and the transmitted wave 
is absent. In the resonance, the reflected wave is totally 
suppressed, the interference pattern in the far field dis- 
appears, and the transmitted wave is clearly seen under 
the layered superconductor. 

Conclusions. — The enhanced transmissivity of modu- 
lated slabs of layered superconductors has been predicted 
for terahertz radiation. It results from the diffraction 
of the incident wave and the resonant excitation of the 
waveguide modes (WGMs). The diffracted wave does 
not decay deep into the superconducting slab, contrary 
to metals, where the enhanced light transmission is due 
to the excitation of evanescent surface waves. A series of 
resonance peaks, associated with excitation of different 
WGMs, is observed in the dependence of the transmit- 
tance T on the incidence angle 9. 
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FIG. 2: (Color online) The transmittance |To| 2 and re- 
flectance |i?o| 2 vs the incidence angle 8 for d — X c — 10 A^f, , 
e = 16, r a6 = 10~ 4 , / = 0.2, fi = 1.1, gc/co = 3. The first 
number in the brackets shows the diffraction order respon- 
sible for the resonance and the second is the number of the 
resonantly excited waveguide mode. 
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FIG. 3: (Color online) The magnetic-field distribution: (a) 
the nonresonance case at 6 = 45°, (b) the resonance diffrac- 
tion for the order n = +1, 8 = 66.1°. The other parameters 
are the same as in Fig. [2] The superconducting slab is mag- 
nified tenfold in the ^-direction. 
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